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We present a quantum theory of nondegenerate phase-sensitive parametric 
^ ' amphfication in a x*^'^^ nonhnear medium. The non-zero response time of the 

^ \ Kerr (x^^"*) nonhnearity determines the quantum-hmited noise figure of x^^^ 

^ \ parametric amphfication, as well as the limit on quadrature squeezing. This 

non-zero response time of the nonhnearity requires coupling of the parametric 
^ ■ process to a molecular-vibration phonon bath, causing the addition of excess 

noise through spontaneous Raman scattering. We present analytical expres- 
sions for the quantum-limited noise figure of frequency non-degenerate and 



Ch ' frequency degenerate x parametric amplifiers operated as phase-sensitive 

! 

^ , amplifiers. We also present results for frequency non-degenerate quadrature 

squeezing. We show that our non-degenerate squeezing theory agrees with 
, , the degenerate squeezing theory of Boivin and Shapiro as degeneracy is ap- 

' proached. We have also included the effect of linear loss on the phase-sensitive 

process, (c) 2008 Optical Society of America 
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1. Introduction 

Fiber-optical parametric amplifiers (FOPAs) are currently the subject of much research for 
use in wavelength conversion^ and efficient broadband amplification.^ They are also candi- 
dates for performing all-optical network functions. Advances in pumping techniques have 
permitted improvements of the noise figure (NF) of FOPAs operated phase-insensitively^'^ 
and the manufacture of high-nonlinearity and microstructure fibers has improved the gain 
slope^'S of FOPAs. 
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In order to explain our experimental noise figure result for a FOPA operated as a phase- 
insensitive amplifier (PIA),^ we have recently published a quantum theory of x^^^ parametric 
amplifiers that takes into account the non-instantaneous nonlinear response of the medium 
and the requisite addition of noise caused by this non-instantaneous nonlinear response. 
The FOPA operates as a PIA when an idler beam is generated wholly in the fiber or when the 
phase difference (twice the phase of the pump less the phases of the Stokes and anti-Stokes 
beams) is approximately or vr and the Stokes and anti-Stokes beams are approximately 
equal in amplitude. Our work with FOPAs operated as PIAs also provides analytical expres- 
sions for the noise figure of x^^^ phase-insensitive parametric amplifiers'^ and wavelength 
converters. This theory shows excellent agreement with experiment.^ In addition, we have 
recently experimentally investigated the noise figure spectrum for PIA and wavelength con- 
verter operation of a FOPA, and shown good agreement to an extended theory that includes 
distributed loss.'^ This inclusion of distributed loss was necessary to model the experi- 
ment,'^ but also provides the necessary quantum theory for predicting the performance of a 
distributed amplifier. 

Phase- sensitive amplifiers (PSA)''^''^ are also of interest because unlike PIAs, they can 
ideally provide amplification without degrading the signal-to-noise ratio (SNR) at the in- 
put.'^ Operation of a FOPA as a PSA occurs when the phase difference at the input is 
approximately 7r/2 and the signal and idler are approximately equal in amplitude. Experi- 
ments with fully frequency degenerate fiber phase-sensitive amplifiers have demonstrated a 
noise figure of 2.0 dB at a gain of 16 dB,'^ a value lower than the standard phase-insensitive 
high-gain 3-dB quantum limit. A noise figure below the standard PIA limit has also been 
reached in a low-gain phase-sensitive amplifier.'^ However, these fully frequency degenerate 
PSA experiments were impaired by guided-acoustic-wave Brillouin scattering (GAWBS)'^ 
requiring pulsed operation'^ or sophisticated techniques for partially suppressing GAWBS.''' 
In order to avoid the GAWBS noise one may obtain phase-sensitive amplification with an 
improved experimental noise figure by use of a frequency nondegenerate PSA. In addition, 
the nondegenerate PSA, unlike its degenerate counterpart, can be used with multiple chan- 
nels of data. A nondegenerate PSA is realized by placing the signal in two distinct frequency 
bands symmetrically around the pump frequency with a separation of several GHz, so that 
GAWBS noise scattered from the pump is not in the frequency bands of the signal. Such 
frequency nondegenerate regime has been demonstrated experimentally, showing good agree- 
ment with theory for the average values of the signal and idler. ^'^ We have also experimentally 
demonstrated an improved bit-error rate by use of a PSA as opposed to a PIA of comparable 
gain.^' However, quantum-limited noise figure measurements have yet to be performed. So 
an analysis of this case is practically useful. Accordingly, we here describe in suitable detail 
a quantum theory of FOPAs that takes into account the non-zero response time of the x^^^ 
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nonlinearity along with the effect of distributed hnear loss. We present analytical expressions 
for the quantum-limited noise figure of CW x^^^ PSAs in the frequency nondegenerate case. 
We also report the limiting value of the NF when degeneracy is approached. In addition, 
we note that this theory is valid for FOPAs used as distributed phase-sensitive amplifiers. 

A frequency nondegenerate parametric amplifier can also operate as a phase-sensitive 
deamplifier (PSD) of two-frequency input signals. This occurs when the phase-difference at 
the input is approximately — 7r/2 and the Stokes and anti-Stokes inputs are approximately 
equal in amplitude. When a PSD is operated with no input signal, such a parametric ampli- 
fier is said to produce "quadrature-squeezed vacuum" (parametric fluorescence of the PSD) 
whose two- frequency homodyne detection exhibits photocurrent variance less than that of the 
vacuum for suitable choice of homodyne phases. Quadrature squeezing has been proposed 
for applications in quantum communications, ^^"^^ improved measurement sensitivity,^^' 
and quantum lithography^^. In the case of FOPAs, previous work by Shapiro and Boivin'^" 
used the dispersionless theory of self-phase modulation developed by Boivin and Kartner ^^'^^ 
that included the non- instantaneous response of the x^^^ nonlinearity to obtain a limit on 
quadrature squeezing in the fully four-degenerate-wave case. In this paper, we present results 
for frequency nondegenerate CW quadrature squeezing for a noninstantaneous nonlinearity 
in the presence of dispersion. We show that optimal squeezing occurs for slightly different 
input conditions than those for optimal classical deamplification. In addition, we show that, 
unlike the dispersionless case, the degree of squeezing reaches a constant value in the long- 
interaction-length limit when the linear phase-mismatch is nonzero and a non-instantaneous 
nonlinear response is present. Our nondegenerate squeezing theory agrees with the previous 
degenerate squeezing results of Boivin and Shapiro'^'' when degeneracy is approached. 

This paper is organized as follows: In Section 2 we discuss the solution of the equations 
describing evolution of the mean values of the pump. Stokes, and anti-Stokes fields. In Section 
3, we present a quantum mechanically consistent theory of the FOPAs. For calculation of 
the noise figure in phase-sensitive operation, we need to obtain only the mean and variance of 
the total photocurrent at the Stokes and anti-Stokes wavelengths. Thus we chose to calculate 
only the output Heisenberg annihilation operators for the Stokes and anti-Stokes frequency- 
pair of interest. We thus obtain the desired results in a simpler way than had we used the 
Wigner or positive-P formalism developed by Drummond and Corney^'^ for propagation of 
the quantum states of pulses in the fiber. In Sections 4 and 5, we apply this theory to obtain 
the noise figure of phase-sensitive amplification and to obtain the degree of nondegenerate 
quadrature squeezing, respectively. We reemphasize the main results and conclude in Section 
6. 
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2. Classical phase-sensitive amplification and deamplification 

We have discussed the x^^^ nonlinear response at length; only a very brief summary is 
presented here. 

The nonlinear refractive index of the Kerr interaction can be written as 

4eongc' 

where Uq is the linear refractive index of the nonlinear material, eo is the permittivity of free 
space, and c is the speed of light in free space. For clarity, we state that 

,2" 



A — Aim 



m 

V2 



(2) 



The x^^K't) nonlinear response is composed of a time-domain delta-function-like electronic 
response (<^ 1 fs) that is constant in the frequency domain over the bandwidths of interest 
and a time-delayed Raman response (~ 50 fs) that varies over frequencies of interest and 
is caused by back action of nonlinear nuclear vibrations on electronic vibrations. Recent 
experimental and theoretical results demonstrate that the nonlinear response function x^^^ (^) 
can be treated as if it were real in the time domain, ^^'^^ yielding a real part that is symmetric 
in the frequency domain with respect to pump detuning and an imaginary part in the 
frequency domain that is anti-symmetric. 

Although a nonlinear response is also present in the polarization orthogonal to that of the 
pump, this cross-polarized nonlinear interaction is ignored because we assume that the pump. 
Stokes, and anti-Stokes fields of interest stay copolarized as their polarization states evolve 
during propagation through the FOPA. Parametric fluorescence and Raman spontaneous 
emission are present in small amounts in the polarization perpendicular to the pump, but 
do not affect the NF of the amplifler. 

We can write A^2(^) in the frequency domain as a sum of electronic and molecular contri- 
butions: 

N^in) =n2c + n2rF{n). (3) 
We also deflne a nonlinear coefficient 7q to be 

27rA^2(^^) 
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where A is the pump wavelength and A^s is the fiber effective area. Thus our 70 is equivalent to 
the nonlinear coupling coefficient 7 used in Agrawal."^^ It is the scaling of A^s with wavelength 
that mainly causes 7^ to be no longer anti-symmetric with 7_q at detunings greater than 
several THz. In what follows, our analytical treatment of the mean fields allows for the more 
general case of asymmetry in the Raman-gain spectrum. However, other results including 
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graphs assume an ant i- symmetric Raman spectrum as this has a minor effect on the quantum 
noise at large detunings. 

We next present solutions to the mean field equations governing a parametric ampli- 
fier. The optical fields are assumed to propagate in a dispersive, polarization-preserving, 
single-transverse-mode fiber under the slowly-varying-envelope approximation. As the in- 
volved waves are quite similar in frequency, to good approximation all fields can be treated 
as if their transverse mode profiles are idententical. Even though the fibers used to construct 
FOPAs typically support two polarization modes and the polarization state of the waves is 
usually elliptical at a given point z in the FOPA, for typical fibers it is still appropriate to 
describe the system with a scalar theory if the detuning is relatively small. This is because 
the input waves are copolarized at the beginning of the amplifier and the fields of interest 
remain essentially copolarized during propagation down the fiber. 

Consider the field 

A{t) = A^^ As exp(ifit) + Aa exp{-int) (5) 

for the total field propagating through a FOPA having a frequency and polarization degen- 
erate pump. Here 1^4^ has units of Watts. The lower frequency field we refer to as the Stokes 
field. As] the higher frequency field referred to as the anti-Stokes field, Aa- The classical 
equation of motion for the total field^^ with the addition of arbitrary frequency dependent 
loss is: 



dA{t) 
dz 



J dT-f{t-T)A*{T)A{T) A{t) - J ^^A{n)exp{-int)dn, (6) 



where a{Q) is the power attenuation coefficient at detuning Q from the pump and A{Q) is the 
Fourier transform of the field. Because the involved waves (Stokes, anti-Stokes, and pump) are 
CW, the usual group- velocity dispersion term does not explicitly appear in Eq. (6). However, 
dispersion is included; its effect is simply to modify the wavevector of each CW component. 
Taking the Fourier transform of Eq. (6) and separating into frequency-shifted components 
that are capable of phase-matching, we obtain the following differential equations for the 
mean fields:^^ 

^ = ^lo\A,U,-^A,, (7) 

^ = ii^o + ln)\Ap\^Aa + t^nAlA:expi-tAkz)-^Aa, (8) 

^ = i{io + l^n)\Ap\^As + t^^nAlAlexp{-iAkz)-'^As. (9) 
d^; ^2 

Here Ak = ka + kg — 2kp is the phase mismatch. Expanding the wavevectors in a Taylor 
series around the pump frequency to second order, one obtains Ak = /?2^^ to second order. 
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where (32 is the group- velocity dispersion coefficient. The attenuation coefficients are aj for 
j = p,a,s at the pump, anti-Stokes, and Stokes wavelenths, respectively. The nonlinear 
coupling coefficients 70, 7n, and 7_n are as defined in the previous section. Eqs. (7-9) are 
valid when the pump remains essentially undepleted by the Stokes and anti-Stokes waves 
and is much stronger than the Stokes and anti-Stokes waves. The solution to Eqs. (7) and 
(9) can be expressed as 

Aa{z,L) = fla{z,L)Aa{z) + Ua{z,L)Al{z), (10) 

As{z,L) = fis{z,L)As{z) + Us{z,L)A:{z), (11) 

where we have explicitly written the solution as a function of both a starting point z for the 
parametric process and an end point L for the fiber. We do this because we will be interested 
not only in the input-output relationships of the electromagnetic fields, i.e. ^^(0, L), but also 
evolution of noise generated at a point z that propagates to the end of the fiber, L. In the 
following subsections, we provide expressions for fij{z,L) and Vj{z,L) for the three main 
cases of interest. 

2. A. Distributed loss solution 

In the most general case, when there are no restrictions on A/c and distributed linear loss 
is present, Eqs. (7-9) can be shown to have a series solution. We here briefiy outline the 
derivation of this solution. Solving for the mean field of the pump, Eq. (7), we obtain 

A^{z) = exp[^7o/p(0)2eff - ^]^p(0), (12) 

where the effective length z^s is defined to be Zcs = [1 ^ exp{—apz)]/ap. Further defining 
the initial pump power in Watts to be /p(0) = |y4p(0)p, and setting the reference phase to 
be that of the pump at the input of the fiber, we substitute the resulting expressions into 
Eqs. (8) and (9). Writing 

Aa = Baexp[i{'jo + -fn)lp{0)zcs - aaz/2], (13) 
A, = Bsexp[i{-fo + -f^n)Ip{0)z^s-asz/2], (14) 

and making a change of variable from z to z^s, one obtains 



d^eff 

di?: 



d2;efr 

where 



^7o/p(0)exp[-/(^efr)]5:, (15) 
-z-f*_M^)exp[f{z,f,)]Ba, (16) 



X .r , * lr^n^ (tts - + 2zAfc) ln(l - apZ^g) 

fiZes) = i[ln + l-Q]Ip{OjZcS ^ • (17) 
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After some algebra and making use of the substitutions 

Fa = 5,exp[/(zeff)], (18) 
F: = i?:exp[-/(zeff)], (19) 

we can obtain the nonhnear coupled equations 

^-[r + -^^^ — ]Fa = (20) 
^ + [r + ^ ]f: = 6^,, (21) 

where the following constants are used for calculating evolution from point z to L: T = 

^bn + 7-f7]^p(^)/2, A = [as/2 - a„/2 + iAk]/2, = i-fnlp{z), and 6 = -z7-q^p(^)- Using 
the expansion 

1 



1 OipZeS n=0 

on the nonlinear term in Eqs. (20) and (21), we find a series solution for Fj and then obtain 
A,. 

The series solution converges in relatively few terms when apZ^s is small, which is the case 
for practical amplifiers. The following solutions are then obtained: 

oo 

fia{z, L) = exp[p{z, L)] J2 (^nL'^s («o = 1, Sg = 0), (23) 

n=0 

oo 

/i,(2;,L) = exp[p{z,L)]Y.SnL:f, (ao = 0,s* = l), (24) 

n=0 

oo 

Uaiz, L) = exp[p(2;, L)] ^ a„L^g (oq = 0, Sq = 1), (25) 

n=0 

oo 

i^siz, L) = exp[p(2;, L)] ^ s„L^fj (oq = 1, Sq = 0), (26) 

n=0 

where 

piz,L) = t['yo + {ln-<n)/'^]WLeE (27) 
-iAk{L - z)/2 - aa{L - z)/4 - a,(L - z)/4 (28) 

and Lcfj = {1 — exp[— ap(L — z)]}/ap. The coefficients a„ and s* are then calculated through 
the following recursion relations: 



(29) 
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2.B. Lossless, Ak ^ solution 

The solution for lossless fiber and Ak 7^ is well known, as are the n and v functions which 
can be expressed as"^^ 



fia{z,L) = exp 



z AA;- [27o + 7o-7-n]/p {L - z 



IK 



X I — sinh[5f(L — z)] + cos\i[g{L — z)] 



(31) 



Hs{z,L) = exp 



z (AA; - [270 + 7-n - Ynllp) {L - z) 



IK 



— sinh[^*(L - z)] + cosh[^*(L - 2)] 1 , 



(32) 



..(.,L) ^ exp I - ^ T "-"^'-^ " s.nh[,(L - .)], 



(33) 



i^s{z,L) = exp 



t {Ak - [270 + 7-n - 7n]^p) - 2) \ n-f7^p(^)' 



9* 



sinh[5f*(L — z)]. 



(34) 



Here Ip{z) = Ip = \Ap{0)\'^ is the pump power in Watts, k = Ak + (7^ + ■ylQ)Ip, and 
g = J —{k/2Y + 7n7*Q-^p is the complex gain coefficient. 



2. C. Lossless, Ak = solution 

We also state the results for the lossless, Ak = case, which is useful in our analysis near 
degeneracy. We have Ak = when the FOPA is pumped at the zero dispersion wavelength 
or if the system is treated as if dispersionless. When the three frequencies are very nearly 
degenerate, we can also make the approximation that Ak = 0. Then the /i and u functions 
become 



^■a{z,L) 

Hs{,z,L) 

Va{z,L) 

Us{z,L) 



exp 
exp 
exp 
exp 



270 


+ 7f^ - 


7*n 


270 


2 

+ l-n 


-7^ 


270 


2 

+ ln- 


7-n 


270 


2 

+ l-n 


-7^ 


2 



Ip{L - z) 


[1 + i-^ni^L - z)Ip], 


(35) 


Ip{L - z) 


[1 + i'y-n{L - z)Ip], 


(36) 


Ip{L - z) 


ijn{L - z)Ap{zf, 


(37) 


Ip{L - z) 


i'j-niL - z)Ap{zf. 


(38) 



Under the normal assumption of an anti-symmetric Raman gain profile (7^ = 'jIq), we 
see the power gain for an input only on the anti-Stokes side, Ganti-stokes = |/^a(0,lv)p = 
1 — 2lm{'~fQ} IpL + |7j7p/pL^, has a four wave mixing gain that is quadratic as a function of 
fiber length and that the Raman loss at the anti-Stokes wavelength is linear in pump power 
and length {—2lm{'~fn}IpL). Similarly, the Raman gain at the Stokes wavelength is linear in 
power and length (2Im{7_Q }/pL). 

2.D. Optimal classical phase-sensitive amplification and deamplification 

We next find the optimal phase-sensitive amplification and phase-sensitive deamplification of 
a mean field consisting of a superposition of Stokes and anti-Stokes fields. We define optimal 
phase-sensitive amplification (deamplification) as the greatest (least) output signal power 
possible for a fixed amount of input signal power. Assuming coherent signal inputs Q having 
powers ICjp and phases exp(z6'j) for j = a, s, the phase-sensitive gain of the PSA is 

^ \A^{L)\' + \ML)\' 
|A,(0)P+|^(0)|2 







ICal 


' + ICI 


2 



M + k.niCaP + (kaP + M\Cs? + [[{l^sK + l^aK)\Ca\\Q e^^'^^'^^ + C.C.} 

\Ca?+\Cs? 

(39) 

By properly choosing the relative power of the Stokes and anti-Stokes inputs and their sum 
phase 9 = 9a + 0s relative to the input pump phase, one achieves maximimum (minimimum) 
phase-sensitive amplification (deamplification). The optimum sum phases 6'psa, opt and 6'psd, opt 
are 

6'psa, opt = - arg[/i^z/; + /i^z/*], (40) 
6'psd, opt = vr - aig[^isK + ^^aK]^ (41) 

for amplification and deamplification, respectively. By setting the sum input power \C,a\^ + 
ICsP = C to be some constant C, the extrema of Eq. (39) can be found to occur when the 
proportion of input anti-Stokes power to the total input power is 

\Ca? =\(l± l/^sP-KP 

where the negative root corresponds to optimum phase- sensitive amplification, the positive 
root to optimum phase-sensitive deamplification. The maximum PSA gain, GpsA is found 
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by insertion of Eqs. (40) and (42) into Eq. (39). The result simplifies to: 



<-^PSA — 



2^2 



2 

(|/i,p-|/i,p)(|z/,p-|z.,p) 



(43) 



In Fig. 1 all of the plots are optimal in the sense that the best total phase and relative 
power of the input fields is chosen. It is clear that the Raman effect is negligible as those 
curves including the Raman effect (circles and squares) are very similar to those neglecting it 
(solid curve and dashed curve). This plot shows also that the definition of the effective length 
is a mathematical one and not a good guide for estimating the gain profile. The dotted and 
dash-dot curves show the gain spectrum of a lossy fiber (0.41 dB/km) 4.44 km in length. 
The other curves are for a lossless fiber of effective length Leg = 3.63 km. Thus distributed 
loss has a greater effect than might be supposed: fibers experience noticeably less gain than 
lossless fibers do when the lossless fiber has a length equal to the effective length of a lossy 
fiber. 

Some of the characteristics of the classical phase-sensitive response can be seen in Fig. 2, 
which is a plot of the phase-sensitive gain vs. fiber length. The primary feature of phase- 
sensitive amplification is that the mean-field gain is relatively insensitive to the relative 
strength of the two input fields (this can be seen by the overlap of the thin solid lines with 
the squares and circles). In addition, typical distributed losses do not significantly impact 
the gain of the fiber, as can be seen by comparison of the squares (lossless fiber) with the 
circles, which represent fiber with loss of 0.25 dB/km. On the other hand, the achievable 
degree of phase-sensitive deamplification is much more sensitive to the relative proportion of 
the input fields which can be seen by comparison of the dash-dotted and dashed lines (equal 
power splitting, i.e., |CaP = ICsP) with the dotted and thick solid line (optimum relative 
proportion). In addition, distributed losses also set a limit on classical deamplification as 
can be seen by comparison of the dashed line (lossy) with the dash-dotted lines (lossless) 
and comparison of the dotted line (lossy) with the thick solid line (lossless). 

3. Input-output quantum mode transformations 

In this section, we discuss the quantum mechanics of the x''^^ parametric amplifier with a 
strong, undepleted, coherent-state pump and derive input-output mode transformations in 
the Heisenberg picture which can be used to calculate the noise figure of the phase-sensitive 
operation of a FOPA and the accompanying quadrature squeezing. Here we also extend our 
previously described quantum theory^"'^^ to include the effects of loss. In order to make 
our treatment consistent with the customary formalism in quantum optics, we rescale our 
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field to a photon-flux field, i.e. the expectation of the field {A) = A = Ajyhu) so that \A\^ 
has units of photons/sec. We also rescale our nonlinear coefficient = '^qTiuj, where % is 
Planck's constant over 2ti. 

We begin with the continuous-time quantum equation of motion for a multi-mode field in 
the presence of a non-instantaneous nonlinearity. This model was presented by Carter and 
Drummond,^^'^^ solved for the case of dispersionless self-phase modulation by Boivin,^^ and 
also derived in detail by Kartner."^^ We have 



dA{t) 
dz 



dr7(t-r)i"^(r)i(r) 



A{t) 



+ im{z,t)A{t)+l{z,t), 



(44) 



wherein 7(t) is the causal response function of the nonlinearity, i.e., the inverse Fourier 
transform of 7^ in Eq. (4). In Eq. (44), rh{z,t) is a Hermitian phase-noise operator 



m{z,t) = HdQ V 2^"'^^"} 1^ ^l^^^y nt _ ^ dMe-'""'}, (45) 
JO Zn 

which describes coupling of the field to a collection of localized, independent, medium oscil- 
lators (optical phonon modes). In addition, l(t) is a noise operator 



lit)= dn^exp{int)vn{t), (46) 

which describes the coupling of the field to a collection of localized, independent, oscillators 
in vacuum state. This coupling is required to preserve the continuous-time commutators 

[A{t),A^it')]=S{t-t'), (47) 

[i(t),i(t')] = 0. (48) 
Note that the time t is in a reference frame traveling at group velocity Vg, i.e., t = 



^stationary frame " ^- The Weighting fuuctiou -y/2Im{7n} must be positive for f2 > so that 
the molecular vibration oscillators absorb energy from the mean fields rather than providing 
energy to the mean fields. The operators dn{z) and d^^z) obey the commutation relation 

[dn{z)Jl{z')] = 6{Q - Q')6iz - z') (49) 

and each phonon mode is in thermal equilibrium: 

{dl{z)dn'{z')) = 6{n - n')6{z - z')nt^ (50) 
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with a mean phonon number rith = [exp(^f2/A;T) — 1]"^. Here h is Planck's constant over 27r, 
k is Boltzmann's constant, and T is the temperature. The operators corresponding to vac- 
uum modes mixing into the Stokes and anti-Stokes frequencies, v±n, obey the commutation 
relations 

[v±n{z), viniz)] = S{±n - ±Q')6{z - z') (51) 

and have no photons in them, i.e., {v±q{z)v±q/{z')) = 0. 

We assume that the total field present at the input of the fiber contains only a single- 
frequency pump, a Stokes field, and a symmetrically placed anti-Stokes field. The operators 
corresponding to these modes are 

A{t) = Ap + As exp(ifit) + Aaexp{-int). (52) 

We do not need to consider other modes in the fiber when the pump is strong (|^pP ^ 
l^sp) and remains undepleted and the Stokes and anti-Stokes frequencies of interest 
are too weak to serve as pumps to other nonlinear processes. When this latter condition 
occurs, fiuctuations from other higher order mixing modes are not coupled in. In addition, the 
vacuum modes at the input do not grow due to parametric fiuorescence to become sufficiently 
strong to serve as pumps for nonlinear processes. This is because the input Stokes and anti- 
Stokes fields would saturate the pump long before this would happen. Thus the symmetric 
pairing of the Stokes and anti-Stokes fields of interest is justified and the pairs do not couple 
through the nonlinear process to other frequencies. 

In order to obtain the coupled- wave equations for the three frequencies of interest, we insert 
Eq. (52) into Eq. (44), take the Fourier transform of the resulting equation, and separate 
into different frequencies. The resulting nonlinear quantum operator equations are: 

i^oAlApAp + ^n]AlAsAp + z[7o + ^^n]AlAaAp 
i[l-n + "yn] AlAsAaexp{iAkz) - ^jThn^jAadl^iz) exp{-i{kp - ka)kz) 
+ ^2lm{^a}AsdQ,{z)exp{-i{kp - ks)kz) + ^j2lm{-fn_^o}Ap[dn~,o+{z) + dn^o.{z)] 
v/a^^^p(^), (53) 

i[lo + ln]AlApAa + i%AlAaAa + i[7o + l2n]AlAsA, 
i^nAlAlexp{-iAkz) - ^j2lm{^n}Apexp[i{kp - ka)z]dn{z) 

Y^a + y/a^Vaiz), (54) 
i[lo + l-n]AlApAs + t^oAlAsAs + ^[7o + 7-2n]il^i. 
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- -^As + y/alvs{z). (55) 

Eqs. (53 - 55) are nonlinear in the quantum operators (except for the last two terms of each 
equation) and are difficult to solve. However, one may linearize these equations around the 
mean values of the operators and obtain solutions accurate to first order in the fiuctuations. 
Using the definition 

Aj = A,+aj, (56) 

where Aj represents a c-number mean value of Aj and dj represents the quantum fiuctuations 
of Aj, and where j = p, a, s, we expand Eqs. (53 - 55) to obtain the following equations for 
the mean values (i.e., those terms containing no fiuctuation operators): 

dAr, 



dz 



dA^ 
dz 



dA 
dz 



i%A;ApAp + i[7o + ^n]AlAsAp + ^[70 + ^_n]AlAaAp 

— — — — 

+i[%n + ^n]A*pAsAa exp{iAkz) - -^Ap, (57) 

^[lo + ln]A*pApAa + i^oAlAaAa + z[7o + l2n\AlAsAa 

— — (y, — 

+i^nAlA: exp{-iAkz) - ^, A^ (58) 

^[70 + l^n]A*pApAs + i^A^AsAs + i[7o + ^^2n\AlAaAs 
+t^.nAlA: exp{-iAkz) - y ^. (59) 



The equations for those terms which contain one fiuctuation operator are 



dd 



^ - i^o[2\Ap\^ap + Alal]+i[^o + ln][\As\''ap + A*Apa^ 



dz 



p pi ' L 'U ' 'S'JLI s| ^p I ^ -^s P s I ^ '■S-' '■p 

+^[7o + l-n][\Aa\^ap + A^Aptta + 

+i[^-n + 7n][^;^da + A^AaUs + A.^dJ] exp{iAkz) 

-^J2lm{'yQ}Aadl^{z) exp{-i{kp - ka)kz) + ^j2hiv{^}AsdQ{z) exp{-i{kp - ks)kz 

a 



ddg 
dz 



dhs 
dz 



+ v2Im{7Q^o}^pMf7-^o+(2;) + d\i^Q-{z)] - ydp + ,Ja;Vp{z), (60) 

^[70 + lQ\[\Ap\'^aa + A^AaCLp + Ap^dJ] + Z7o[2|Aapda + A^d^] 

+«[7o + l2u][\As\'^da + AlAaCLs + A,A„dJ] + iin[Ala\ + 2ApAlap] exp{-iAkz) 

j ^ cx 

+ y2Im{7Q}Apexp[i(A;p - ka)z\dn{z) - yd^ + ./a2va{z), (61) 

^[7o + l-u][\Ap\^ds + AlAsdp + ApAsdl] + z7o[2|A,|2d, + A^dl] 

+i[lo + l~2n][\Aa\^ds + AlA.da + AaA,d\] + i^_n[Ald\ + 2ApAldp] exp{-iAkz) 

a. 



2Im{7n}ipexp[z(/cp - ks)z]dl^{z) - yd^ + y/a^Vs{z). (62) 
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When the amphfier is operating in the unsaturated regime, the pump is depleted to a negh- 
gible degree. Thus one may neglect the four wave- mixing terms in Eq. (57). We next make 
the strong pump approximation, i.e., \Ap\'^ ^ l^aP? l^sP ; which is valid at the input of the 
fiber under typically used operating conditions and remains valid when there is essentially 
no pump depletion. Under these conditions, we note that the cross-phase modulation of the 
pump wave due to the Stokes and anti-Stokes fields is negligible. We thus obtain the final 
form of the mean-field equations: 

^7ol4l'4- (63) 

— — — — Qi — 

i (70 + 7f^) \Ap\^Aa + nnAlAl exp {-lAkz) - -jA^, (64) 

— — — — (y — 

i (70 + 1-n) WAs + ij-nA-lAl exp {-iAkz) - y A,. (65) 

The solution of these mean-field equations [Eqs. (7-9)] were examined in the previous sec- 
tion. Examining the fluctuation equations, we note that the first-order fluctuation equations 
[Eqs. (60-62)] are linear in their quantum operators. The higher-order fluctuation terms are 
nonlinear in the quantum-fluctuation operators, but their contribution is negligible because 
they cannot contain terms proportional to A^. We thus neglect the higher-order fluctuation 
equations which are not reported here. 

Similarly, those fluctuation terms in Eqs. (60-62) that do not contain two pump factors 
can be neglected, as they are much weaker than the other terms. We remark here that the 
Langevin-noise terms due to the Raman effect in Eqs. (61-62) will not be neglected. This is 
justified because comparison of the amplitude of the Langevin terms, ^J21m{^}Ap, in Eqs. 
(61-62) with the the amplitude of those four- wave- mixing terms in Eqs. (61-62) that contain 
only one pump term, for example j-nApAl^ap shows that the Langevin term is much greater 
than the one-pump-field four-wave-mixing terms. For example, consider the symmetrized 
magnitude of the two terms: 



2Im{7_o}|Ap| 


2 


[{dldn + dndl)]/2 _ 


Im{7_Q}(nth + 1) 


47^n 




'\Ap 


\'^[{alap + apal)]/2 \ 


2|7-n 





(66) 



For an anti-Stokes amplitude of up to 10 mW (a flux of Ua = 7.8 x 10^^ photons/s ) this 
ratio is greater than 25 for pure silica and even better for materials with stronger Raman 
gains (such as dispersion-shifted and highly nonlinear fibers). This means that the Langevin 
noise term should be included. 

One last simplification can be made to the pump equation, Eq. (60). We simply note that 
all of the remaining noise terms in Eq. (60) affect only the pump field and that under the 
undepleted, strong-pump approximation, the Stokes and anti-Stokes fields do not interact 
with the pump fluctuations to first degree. This is the case if one assumes that the pump 



dAp 
dz 

dAg 

dz 
dA 
dz 
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is, or nearly is, in a coherent state. Thus, as we are only interested in the noise introduced 
during the amplification process to the Stokes and anti-Stokes fields, we neglect the pump 
fiuctuations and replace Eq. (60) with 

dd^ 



0. 



(67) 



This allows one to treat the pump field fully classically. To summarize the steps taken, the 
undepleted, strong, coherent-state pump approximation yields the following equations for 
the first-order fiuctuations: 
doa 



dos 

dz 



an 



-fY'2Im{7n}Apexp[z(A;p - ka)z\d^{z) - —CLa + s/a'aVa{z), 
i[7o + 7_Q]|Appa^ + i7_QApa^exp(-iA/c2;) 



-Y'2Im{7n}y4pexp[z(/cp - ks)z\dl^{z) - —d^ + ^/alvs{z). 



(68) 



(69) 



We note that Eqs. (67-69) can be combined with Eqs. (63-65) by using Eq. (56) to yield: 



cL4p 
dz 
dAa 



"a 



(70) 



(71) 



(72) 

This final form of the equations allows one to evolve the mean fields of all three frequencies 
and the fiuctuations of the Stokes and anti-Stokes wavelengths simultaneously. 
The solution of Eqs. (71) and (72) is 



i (7o + m) \Ap\^Aa + z-fnA^Alexp {-iAkz) - -jA^ 
y^2Im{7n } exp [i ( kp - ka)z]dn{z) + y/a^Va{z), 

^ = i i^o + ^-n) lApl"^ As + i^-nAlAl exp {-iAkz)-^ As 
2lm{^^}Apexp[i{kp - ks)z]dl^{z) + y/a^Vs{z). 



Aa{L) 



/i,(0,L)i„(0) + z/,(0,L)il(0) + 

L 



As{L) 



2Im{7Q} dzAp{z)exp[i{kp-ka)z][ij,a{z,L)-Ua{z,L)]dQ{z), 

dz [^/a2fla{z, L)va{z) + ^/olvaiz, L)vl{z)] (73) 
/i,(0,L)i,(0) + i/,(0,L)il(0) + 
2Im{7n} / dz Ap{z)exp[i{kp- ks)z\[-^s{z,L) + Vs{z,L)]d\^{z). 

•J 

(74) 



+ dz [y/cTsHsiz, L)vs{z) + ^JoTaVs^z, L)v\{z)\ 
Jo 
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In the notation used in this paper, the functions fij{z,L) and Vj{L^z) denote evolution 
from a point z in the fiber to the end of the fiber (L) where the intensity and phase of the 
pump at point z must be used. 

In this section, we have presented a thorough derivation of the input-output mode trans- 
formations that govern a x''^^ parametric amphfier. In the following two sections, we use 
these input-output mode transformations to obtain the noise figure of x*-^"* phase-sensitive 
parametric amplifiers and the squeezing parameter for quadrature squeezing. 

4. Noise Figure of phase-sensitive amplification 

In this section we discuss the noise figure of phase-sensitive amplification, which is defined 

as 

We assume for our treatment here that the difference between the Stokes and anti-Stokes 
frequencies exceeds the bandwidth of the detector. Thus beat frequencies of these two waves 
will not be detected and can be neglected. The mean photon flux at the Stokes and anti- 
Stokes wavelength is {fij) = {Aj{0)Aj{0)) = and is assumed to be in a coherent state. 
In what follows, we neglect the small frequency difference between oOa and Ug. Thus the input 
SNR can be written as 

- {AhD + {AnD - ICaP + ICP " ■ ^ ' 

Calculating the output SNR in a similar way and plugging into Eq. (75), the NF for phase- 
sensitive amplification can be expressed as 

= (KTW^ ' ^''^ 

where the mean output power at each wavelength. Pa and Pg, is 

Pa = |/ianCa|' + kanC|' + (CaC/ia^^:+C.C.), (78) 
Ps = |/i.nC|' + |i^snCa|' + (CaC/isi^:+C.C.). (79) 

In Eq. (77), we have expressed the variance of the output photocurrent as the sum of a phase- 
insensitive portion, (Ara|j), and a phase-sensitive portion, (Aripg), which are calculated to 
be 

(A^i) = PaBa + PsBs, (80) 
(An|g) = 2Q*B^ + 2QB2, (81) 
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where the quantities 

Bj = |/i,f + |z/,f + (2nth + l)|r,f + |c,-ip + l92|', (j = a, s) (82) 

Q = ifiaCa + l^aOif^sCs + iysC), (83) 

Bi = Cxi + r^{ntY, + 1) + jiai's, (84) 

B2 = ^2 + r>th + /^>a, (85) 

have noise terms defined as follows: 

|r,|2 = 2lm{^n} r dz\Ap{z)\^\fXa{z,L)-i^aiz,L)\\ (86) 
Jo 

|r,|2 = -2Im{7_n} / dz\Ap{z)\^\fisiz, L) - Us{z, L)\\ (87) 

Jo 

rL 

|Ca(5)lP = dzaa(s)\fJ'ais){z,L)\'^, (88) 

/•L 

|Ca(.)2r = / dzas(a)Wais)iz,L)\'^, (89) 







= 2Im{7n} / dz Ap{z) exp{-iAkz) 

x[^a{z,L) - iya{z,L)][-fi,{z,L) + Us{z,L)], (90) 

Cxl = "a / dz fia{z,L)Us{z,L), (91) 

Jo 

c*2 = as dziya{z,Lyjj,s{z,L)*. (92) 

JO 

In the above expressions, |ra(s)p represents the integrated amplified noise at the anti- 
Stokes (Stokes) wavelength seeded by thermally populated optical phonon modes that are 
coupled in by the Raman process. The terms |ca(5)ip represent integrated amplified noise at 
the anti-Stokes (Stokes) wavelength seeded by vacuum noise mixed in through distributed 
loss at the anti-Stokes (Stokes) wavelength, while the terms |ca(s)2p represent amplified noise 
at the anti-Stokes (Stokes) wavelength seeded by vacuum noise mixed in through distributed 
loss at the Stokes (anti-Stokes) wavelength. 

In addition, the phase-sensitive terms fia^^s and represent amplified phase-sensitive 
noise seeded by the vacuum noise at the anti-Stokes and Stokes wavelengths. The quantity 
represents amplified phase-sensitive noise seeded by the thermal-phonon fields due the Raman 
effect, and Cxi and Cx2 represent the amplified phase- sensitive noise seeded by the vacuum 
noise due to distributed linear losses. Phase-sensitive noise is present when the photocurrent 
variance with both Stokes and anti-Stokes waves impinging on a detector is different from 
the sum of the individual noise variances of the Stokes and anti-Stokes frequencies. 

4- A. Degenerate Limit 

By taking the limiting value of the NF as f2 ^ 0, we find the NF performance of a fully degen- 
erate FOPA. We find this limiting value of the NF by expanding the anti-symmetric imagi- 
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nary part of 7^ in a Taylor seris and expanding the exponential in nth = {exp[^f2/(/cT)] — 1}^^ 
before allowing 1^ ^ 0. We also use the fact that in this limit, Ak also approaches and the 
optimum power splitting ratio approaches 0.5. This NF limit is: 



4fcT7- (0) 

NFpsA, n^o = 1 H 



1 _ g^NL 



(93) 



where 0nl = 7o|^pP-^ is the nonlinear phase shift and 7^(0) is the slope of the imaginary 
part of 7q as f2 — > 0. We observe that the PSA noise figure for ^ increases to a 
maximum of slightly more than dB and then decreases again and approaches dB in the 
high-gain limit. This unusual behaviour of decreasing noise figure vs. nonlinear phase shift is 
due to the relative scaling of the Raman and FWM processes when A/c = and fact that the 
mean Raman gain of the Stokes and anti-Stokes frequencies vanishes. The total Raman noise 
scales linearly whereas the two-frequency signal undergoes quadratic gain. Consequently, the 
Raman noise being created at the Stokes and anti-Stokes frequencies does not grow. 

4.B. Results 

In Fig. 3, we plot the noise figure vs. PSA gain for several values of detuning for a typical 
highly-nonlinear fiber with a loss coefficient corresponding to 0.75 dB/km. The plots show 
that for detunings achievable by use of electrooptic elements (40 GHz detuning, phase- 
matched, solid curve), the results are almost exactly the same as would be achieved in the 
limit of zero detuning ( dotted). For these simulations, we use realistic values for 7(0) and for 
the distributed loss. We have additionally assumed that the highly-nonlinear fiber has the 
same ratio of Im{7(fi)} to Re{7(f2)} as the standard dispersion-shifted fiber, i.e., the two 
have the same germanium content. Unsurprisingly, we see by comparing the thick curves 
(lossy) to thin curves (lossless), that the distributed loss increases the noise figure. We aslo 
see that as the Raman-gain coefficient decreases, the noise figure improves. 

Interestingly, unlike the PIA case. Fig. 3 shows that the PSA noise figure is greater than 
dB as the gain approaches dB. This occurs because the Raman gain and loss processes 
dominate in the early parts of the amplifier (Raman gain and loss are linear in the early 
parts of the amplifier while the four- wave- mixing gain is quadratic), adding noise to both fre- 
quencies, while the mean field undergoes no net gain due to the Raman loss at one frequency 
and the Raman gain at the other. 

In Fig. 4, we show the gain and noise-figure spectrum for a 4 km fiber with fiber parameters 
as described in the caption. This plot shows that the increasing Raman-gain coefficient with 
detuning causes an increasing noise figure. These results show that for realistic optical fibers, 
a FOPA operated as a PSA can achieve a noise figure below 1 dB for detunings up to 1 THz. 
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5. Nondegenerate quadrature squeezing 

When no light is injected into the FOPA, the quantum correlations between the Stokes 
and anti-Stokes modes imply the presence of quadrature squeezing, which is measurable by 
homodyne detection with a two- frequency local oscillator (LO). 

The difference current of the homodyne detector may be written as 

i = biqa + b% + }l.c., (94) 

where H.c. stands for the Hermitian conjugate of the first two terms; Qa and are the 
annihilation operators corresponding to the anti-Stokes and Stokes components of the LO 
beams, which are in a coherent state each with photon flux |q!lojP, relative intensity 



yj = 







2 




P + l 


OiLO,s\ 


2 



(95) 



and phase 6j for j = a, s. 

The squeezing parameter is defined as the ratio of the photocurrent variance with the 
pump on to the photocurrent variance with vacuum input to the homodyne detector (i.e., 
pump off). 

5. A. Lossless fiber with Raman effect 

In this paper, we concentrate on squezzing results for a lossless FOPA. This is because 
increased nonlinear drive will, theoretically at least, overpower linear loss, leading to no hard 
limit on the achievable squeezing. However, as the Raman effect at each z scales with pump 
intensity as does the four-wave-mixing process, the lossless case illustrates a fundamental 
limit on the achievable squeezing. Using Eqs. (94) and (95), we obtain, after some simple 
algebra, for the lossless Raman-active case: 

S = ^^^'^ 

= [1 + 2(|z/„p + |r,| Vh)]ya + {1 + 2[\us\^ + \rsWl + nt^)]}y, 

-h2{[/i,i^a(l + '^th) - fJ'a'^snth] exp[-i{9a + 9s)] + c.c.}yRJ~yl- (96) 

In order to produce the best squeezing, one must choose the best phase and relative intensity. 
Once again, these two choices are independent. In order to choose the LO phases, we note 
that the third term in Eq. (96) has a negative sign when 

da + 9s='K + arg[/i^I/a(l + rith) - /ial^snth]- (97) 

Since the total LO power is conserved {ua + ?/s = 1), we may use this fact to eliminate Us in 
Eq. (96), which then becomes quadratic in Maximizing the magnitude of the third term 
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in Eq. (96) as a function of Ua yields the additional condition for maximal squeezing, 



'1 + nth) - kapnth 



(98) 



Thus the power splitting of the LOs for maximal squeezing is slightly different from 50% 
and also slightly different from that for maximal classical deamplification. 

Use of this optimal two-frequency LO yields the following optimal squeezing result 



5, 



opt 



1 + + \^s? + \ra?ntb + |r^P(nth + 1) 



(99) 



In order to make a connection with previous work, we show that in the limit of degenerate 
operation we reach the same result as obtained by Shapiro. ^'^ By placing the pump at the zero- 
dispersion wavelength of the fiber, our solutions are then identical to those in a dispersionless 
fiber (assuming the higher-order terms in an expansion of (3 are negligible). We thus use 
the expressions in Eqs. (35) to (38) for yUj and Vj. In order to evaluate the limit as the 
detuning approaches zero, we expand the anti-symmetric imaginary part of 7^ as an odd 
power series around 1] = and take the Taylor series expansion of the exponential in nth = 
l/[exp(/;,f2/(A;T)) — 1]. Then taking the limit as ^ 0, the squeezing approaches the limit 
derived by Shapiro^*^ for the fully degenerate case. This limit is 



SoA^ ^ 0) 



^7o 



1 + 



5'NL + 



2fcT7:(0) 
^7o 



1/2 



(100) 



where 0nl = 7o|^pP-^ is the nonlinear phase shift and 7^(0) is the slope of the imaginary 
part of 7q as r2 — s> 0. 

In Fig 5, the main features of this squeezing theory are illustrated for a lossless fiber. 
First, it is clear that the Raman effect degrades the achievable amount of squeezing. In 
addition, it can be seen by comparing the dashed lines (phase- matched) , dash-dotted lines 
(partially phase-matched) and other lines (A/c = 0) that when the Raman effect is included, 
phase-matching leads to worse squeezing instead of improved squeezing predicted by an 
instantaneous nonlinearity model. Finally, by comparing the Raman-included dashed and 
dash-dotted lines, we see that when A/c 7^ 0, precise unequal balancing of the relative power 
splitting of the two LO frequencies is required. When this optimum splitting is achieved, the 
squeezing can be seen to approach a constant value. However, when Ak = 0, the possible 
amount of squeezing is not bounded, and squeezing asymtotically scales as l/f/'NL, which is 
explained by the quadratic scaling of the four-wave-mixing process and the linear scaling of 
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the Raman process. The hard hmit of constant squeezing with increasing nonhnear drive in 
the phase-matched {Ak 7^ and k = 0) case occurs when the excess noise due to the Raman 
effect balances the strength of the squeezing process. 

6. Conclusion 

In conclusion, we have presented a quantum theory of parametric amphfication in a x*^^^ 
nonlinear medium that includes the noninstantaneous response of the nonlinearity and the 
effect of distributed linear loss. We have applied this theory to nondegenerate phase-sensitive 
amplification and deamplification and have found the input conditions for optimal ampli- 
fication and deamplification. We have also found the input conditions for operation at the 
minimum noise figure, which for detunings < 1 THz is predicted to be < 0.4 dB in the 
high-gain limit for FOPAs made from typical dispersion-shifted fibers. We anticipate that 
nondegenerate phase-sensitive amplifiers will produce record noise-figure performance, as 
they allow circumvention of GAWBS noise that is present in the degenerate case. We have 
also presented a theory of non-degenerate squeezing and found the optimal continuous-wave 
local oscillator for a lossless FOPA with non-instantaneous nonlinear response. Our results 
agree with the limit previously found by Shapiro'^" as degeneracy is reached. Away from de- 
generacy and with a nonzero linear phase-mismatch, we have shown that optimal squeezing 
in a dispersive fiber when Ak 7^ reaches a constant limit, unlike the 1/0nl scaling that 
occurs when the linear phase-mismatch vanishes. 
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Fig. 1. Gain spectra vs. detuning for optimized PSAs made from (a) lossless 
fiber of length L = 3.63 km and Im{7n} = (solid curve), (b) same as (a), 
but Im{7Q} calculated for DSF as explained in the text (circles), (c) 4.44 km 
of DSF with = 3.63 km for = = = 0.41 dB/km and Im{7f7} 
calculated for DSF as explained in the text (dotted curve). Gain spectra vs. 
detuning for optimized PSDs made from (d) lossless fiber of length L = 3.63 
km and Im{7n} = (dashed curve), (b) same as (a), but Im{7Q} calculated for 
DSF as explained in the text (squares), (c) 4.44 km of DSF with L^s = 3.63 km 
for = 0.41 dB/km and Im{7Q} calculated for DSF as explained 

in the text (dash-dotted curve). Input pump power is 0.33 Watts, Aq = 1551.16 
nm, pump wavelength is 1551.5 nm, and the dispersion slope is 57 ps/(nm^ 
km). 
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Fig. 2. Gain vs. fiber length for PSA made from DSF for (a) phase-sensitive 
deamphfication with optimum power sphtting in lossless fiber = = 
ttp = 0) (thick solid line), (b) phase-sensitive deamphfication with optimum 
power splitting in a lossy fiber = as = ap = 0.25 dB/km) (dotted line), 
(c) phase-sensitive deamphfication in a lossless fiber with \(a\'^ = ICP (dash- 
dotted line), (d) phase-sensitive deamphfication in a lossy fiber with = 
ICsP (dashed line), (e) phase- sensitive amplification in a lossless fiber with 
I Cap = ICsP (squares) and optimum input power splitting (solid curve under 
the squares), and (f) phase-sensitive amplification in a lossy fiber with |Cap = 
iCsp (circles) and optimum input power splitting (solid curve under the circles). 
Input pump power is 4 Watts, pump-signal detuning is 1 THz, and phase- 
matching is achieved at the input [AA; = — 2Re{7f7}Jp(0)]. 
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Fig. 3. PSA noise figure vs. gain for various detunings for a highly-nonlinear 
fiber. For thick lines, fiber attenuation is 0.75 dB/km at pump, Stokes, and 
anti-Stokes wavelengths; for thin lines, fiber is lossless. VL/2ti = 13.8 THz, 
dashed curves; Q/27T = 1.38 THz, dash-dotted curves; Q/2tc = 40 GHz, solid 
curves; and f2/27r = Hz, dotted curves. Except for dotted curves, phase 
matching at the input [AA; = —2Re{jQ}Ip{0)] is achieved. For dotted curves, 
Ak = 0. The anti-Stokes/Stokes relative phase and power splitting at the input 
is for optimal classical gain. Initial pump power is 340 mW, 7(0) = 9 x 10~^ 
W~^m~^, peak imaginary part of 7^ is 3.5 x 10^^ W~^m~^. Fiber length is 1 
km. 
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Fig. 4. PSA gain and noise-figure spectrum vs. detuning. Anti-Stokes/Stokes 
relative phase and input-power splitting is for optimal classical gain. Initial 
pump power is 300 mW, 70 = 2 x 10~^ W~^m~^, peak imaginary part of 7^ 
is 0.75 X 10~^ W~^m~^. Attenuation is 0.41 dB/km at the pump, Stokes, and 
anti-Stokes wavelengths. Fiber length is 4 km. Aq = 1551.15 nm, Ap = 1555.5 
nm, and the dispersion slope is 57 ps/nm/km^ 
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Fig. 5. Squeezing vs. nonlinear phase shift in a lossless fiber. Lower lines, 
without Raman effect. Upper lines, with Raman effect. Dashed lines signify 
that phase-matching is achieved at the input [Ak = — 2Re{7n}|ApP]. The 
thick dashed line is for optimal LO power splitting, thin dashed lines for equal 
LO power splitting. Dash-dotted lines are for Ak = — (2/3)Re{7n}|App. Thick 
dash-dotted line is for optimal LO power splitting, thin dashed lines for equal 
LO power splitting. In all other lines, Ak = 0. Raman effect neglected, dotted 
line; Raman effect included and equal LO power splitting, marked with x; 
Raman effect included and optimal LO power splitting, circles; CW limit of 
Eq. (100) with 'jl{0)kT/h'jQ = 0.026, thin solid line. Pump-signal detuning is 
40 GHz. 
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